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LIST OF SYMBOLS 


Svmbol 

fiigUsh: 

r«finition 

Units 

A* 

C<dleetw efficiency parameter 

— 


Gross area of a single collector 

ft* 


Gross area of a single collector including effective area 
of related external manifold 

ft* 


Area of manifold tube (or pipe) betwewi a<tjacent 
collector parts, rrDL 

ft* 

B« 

Colleotor efficiency parameter 

BTU 

ft*.hr.® 


Mean value of flow media ^ecific boat 

BTU 

lbm.®P 

D 

Manifold tube (or pipe) outside diametw, not including 
insulation 

ft 


Total insolation in plane of (selector 

ft^.hr 

*ar 

Incident angle modifier parameter obtained from 
ASHRA E 93-77 test method 


L 

Nominal length of manifold between adjacent collector 
inlet (or outlet) ports 

ft 

ril 

Rate of flow of heat transfer fluid 

IbmA^r 

n 

Number of eoUectors placed in parallel in collector array 
bank 

— 

Q 

Heat loss (or gain) 

BTU/hr 

R 

Overall effective thermal resistance of manifold 
insulation covering 

ft*.hr.® 

t 

Temperature 



*A and B denote the firat order ooUeetor efficiency parametm» respectively intercept and 
negative of slope, as determined from an ASHRAE 93*77 teat. 
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U8T OF SYMBOLS (eont'd) 


Symbol Doftnition 

OrMkt 

a Solar absorptivity of eoUeetor abaorbar plata 

1 } Bffioianoy of solar ooUaetor basad on gross araa (not 

* Ineluding ralatad araa of axtarnal manifold) 

T Solar transmissivity of eollaotor glaaing 


Subsoriptst** 

a 

B 

C 

BM 

g 

I 

j 

L 

M 

O 

u 


Rafaranead to it conditions 

Before mixing at an outlet manifold tarmiyius 

Collaetor 

External Manifold 

Rafaranead to gross araa 

Inlet conditions 

Index (e^M J-I refers to first collector in array) 
Losses 

Manifold and mixed 
Outlet conditions 

Useful (e^M is the ireful enerffy output) 


Units 


**Bxample of subseriptsi tj^^j denotes the manifold (M) inlet (1) temperature (t) for eoUacfor 
number one (Jsl). the manifold (M) outlet (O) tampwatura (t) for mixed (M) 

flow just downstream of the outlet of coUeetcM' number two (J=2). 
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Figure No« Title Pije 

1 Typical Flat Plate Solar Thermal Colleotor With External S 

Manifold Required 

2 Plat Plate Sclar Thermal Colleetor With Internal (integral) 3 

Manifold/Header 

3 Typical Collector Array Bank With External Mawifold 6 

4 Thermal Performance Curves With and Without External 13 

Manifolds — Grumman Collector, Model 332A 

5 Thermal Performance Curves With and Without External 14 

Manifolds -- G^«*al Electric Collector, Model TC 100 


1.0 PURPOSE 

The purpose of this dooument is to present the results of en enelytieal determination 
of the effect of heat loss from solar thermal eoUeotor arrap manlfol<EnR on array 
thermal performance. The analysis oomputer program permits a direot eomparison 
of array thermal performanoe with and without external manifol<Hng. 

2.0 SUMMARY 

A simple FORTRAN oomputer program was written to permit the computation of 
the thermal performance of solar thermal collector arrays with and without external 
manif iids. Computations were carried out for arrays constructed from two example 
solar thermal collectors (Grumman Model 332A and General Electric Model TClOO) 
that are to be used in the DOE Solar in Federal Buildings Program. The results of 
the analysis are presented graphically for typical ex;;emal manifold siaes and 
thermal insulations and are compared directly with the collector alone thwmal 
ptfformance (i.e., the ASHRAE 93-77 thermal performance results based on 
collector alone gross area). 

The results indicate significant degradation in collector array thermal performance 
for manifold insulation having thermal resistance hr ^ t * .” F^ values of 2.0 or less. 

Values of R near 7.0 and greater denote little heat lo£%d performance d^ada- 
tion. Note that the determination of an effective coUector area that includes any 
external manifolding and spacing between adjacent collectors is critical to adequate 
dettf mination of array performance. The effect of the assumed effective area on 
array performance is shown graphically for the two cases computed along with the 
effect of the various insulation R values for the external manifold. 

3.0 INTRODUCTION 

Solar thermal collectors have been designed with a variety of inlet and outlet 
configurations. This is especially true in the ease of flat plate and evacuated tube 
designs and has resulted in a variety of manifolding configurations, manifold pipe 
size, matwials and insulations, varying quality of installation workmanship, etc. 

Ideally, the system designer should minimize the amount of external manifolding 
required. The two most important reasons for this are as follows: 1) to reduce the 
physical space required for a collector array; and 2) to minimse manifold heat 
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Uw6S. It Is obvious thst the amount of useful energy ooUeeted per unit area vrfU be 
diminished If tlK^ eoUeetor array wastes spaoe unduly and If external manlfoldB are 
poorly Insulated. Thus, oolleetors using internal (Integral) manifolds are generally 
iweferred when designing large arrays of sudt oolleetors. 

Althoig^ there are many other advantages, as w<dl as some disadvantages, to the 
use of oolleetors having internal manifold eonfigurations, the intmit of thto report is 
to provide a general oomparison of the area and heat loss related effeets of 
externally versus internally manifolded oolleetors. Seetion 4.0 provides an overall 
discussion of the study, inoluding the baaie assumptions, analytieal model, and 
oomputational parameters. Seetion 5.0 presents two example input data sets (two 
oolleetors) and Seetion 6.0 presents the results of the study. General eonOiusions 
from the study are presented in Seetion 7.0. A list of flat plate solar thermal 
o<dleetors that may qualify as having internal (int^al) manifolds is presented in 
Appendix A. Appendix B presents both a sample eomputer input and output and the 
program listing used in the present study. 

4.0 DISCUSSION 

4.1 Intwnal/External Manifolds-Definition. Except for those solar oolleetors that use a 
single seipentine flow path, most flat plate solar thermal eolleetors have a distinct 
internal header /riser eonfiguration (Figure 1). The eolleetor of Figure 1, as 
configured, requires a significant external manifold for it to be eonneeted in parallel 
to another eolleetor, or eolleetors, in an array bank. Thus, in this report the 
eolleetor of Figure 1 will be treated as havir^ an external manifold in spite of its 
header/riser configuration. 

The eolleetor examples shown in Figure 2 may be considered to fit the description of 
internally (int^aUy) manifolded eolleetors if designed as such. Obviously, such 
eolleetors will function adequately only if the headers and eonneotors must be 
adequately sized to maintain balanced flow rates in all of the risers. 

4*2 Scope of the Present Study. The study described in this report provides a simple, 
accurate comparison of the instantaneous thermal effieieneies for real and idealised 
flat plate solar thermal eolleetors with and without internal (integral) manlfoldB. 
The following parameters were varied to evaluate their effect on eolleetor array 
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Of flET MANIFOLD 



INLET MANIFOLD 
(Typical) 


Figure 1. Typical Flat Place Solar Themal Collector 
With External Manifold Required 



Figure 2. Flat Plate Solar Thermal Collector With 
Internal (Integral) Manifold/Header 
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effioieneyt 


1) Manifold Insulation Effective Thermal Resistance, R s o,7S*, 2 and 

’“TTu — 


2) Ku ocr of Parallel Collectors in Array Bank, n - 8. 

3) Average Ambient Temperature, t s 40 ^f. 

fli 

4) Manifold Inlet Tempwature to Array Bank, t^^j » 40.1 °F, 150 ^F, and 220 
®F. 


4.3 Underlying Assymptims of Analysis. A set of assumptions was established to 

provide a simple, but realistically accurate analysis. It is deemed that no single 

assumption seriously limits the results of the current analysis: 

1) All collectors in an array bank shall be identical and shall u^e a liquid as the 
flow media. 

2) The liquid flow rate through each collector shall be identical and equal to ih . 

3) External manifolds— when used— shall be of a constant diameter tube (or pipe) 
and be uniformly insulated. 

4) Manifold heat losses (Qyy|) shall be deter riiined from the !r.put valuer of the 
effective insulation resistance (R) and the difference fn the calculated mean 
fluid temperature, t^, between collector inlet (or outlet) ports and the 
ambient temperature, t_, i.e., 

Cl 

^LM “ I? 

and where A^ is the total tube (or pipe) area between collector ports on the 
same manifold. 

5) The external manifolds receive no heat inout from exposure to the sunlight. 




An R value of 0.75 


ft^.hr®F 

Sir 


approximates that of an uninsulated pipe. 
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6) The number of eolleotors, n, in an array bank shall be 10 or Iom for 
eomputation purposes. 

7) The groes area of the ooUeetor array bank shall be defined as n times the 
effeottve gross area of a single collector. Two separate eolleetor areas shall 
be iigmti a) - the ^ross o^ a single eoUeotor as defined by A8HRAE 
93-77, and b) - the e/iective gross area of a ooUeotor including its 
externally conflgured manifold. These two areas, A^ and are further 
defined in Figure 3. 


4.4 Required input Data. The following input data is required to perform the an^ysis 
for fJetcrmtning the effect of manifold heat lones and/or effective array size on 
eolleetor array bank performance} 

1) Collector efficiency parameters based on gross area (A^ and first order curve 
fit In accordance with A8HRA B 93-77 test results (A, B, and 


»7g = AK^^+ B (t,-g/It 


( 2 ) 


2 ) 

3) 

4) 
9 ) 
6 ) 

7) 

8 ) 
9) 

10 ) 


where B is defined in English units, K^ir^ incident angle modifier, t| 
is the collector inlet (fluid) temperature, and t^ and are defined below. 


Ambient temperature, t (^F) 

cl 

Manifold inlet temperature, tj^j (^F) 


BTU 

Total insolation measured in the plane of the collector, L ( — ^) 

2 ' hr.ft^ 

Collector gross area, A^ (ft^) 

Collector effective gross area including external manifold, A 


(ft*) 


Flow rate through each collector, (Ibm/hr) 

Mean value of flow media specific heat, c^ (BTU/lbm.^F) 
Number of collectors in array bank, n 


Effective external area of manifold tube (pipe) betwem ai^aeant inlet (or 
outlet) ports, A^ (ft^) 
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11 ) 


Bfk>etiva tlwrmal rMlttanoe of the manIfoM Imulatlon, R (hr.ft^^P/BTU). 
Noto thtt ttM arott Ap, mod in aquation (1) mutt eorraspond to ttw propof 
area at dttignated by the pipe ineulation manufaeturer. In this report the 
manIfoM pipe (tube) diameter it UMd in ealeulating the heat lorn area. 

4.5 Governing Bouationt. For timplteity in the current analytit, only linear equatione 
and mean valuee of individual parameters have been utilised. Thus* the analysis 
proeeedt with the following equations and in the general order presented. 


4.5.1 Inlet Temperature for Inlet Manifold and Individual Collcetorsi Given the tempera* 
ture, tj^p of the flow at the start of the inlet manifold^ the heat loat from each 
length of inlet manifold mint be computed. With the assumption that the 
temperature change of the fluid in the manifold is linearly related to the manifold 
arest ambient temperaturet and insulation thermal resistance, then the temperature 
at the terminus of a manifold length may be computed as follows: 


. _ Wl * ~tP <; 

Equation (3) is derived irom an energy balance on a length of manifold, i.e., 

A 

^LMIJ “ IT ' 'a^ 

end 

with the mean manifold fluid tempwature defined as 


By defining the temperature at the terminus of each inlet manifold 

section may be determined by repeated application of Equation (3). Further, the 
inlet temperature of each collector, is assumed to be equivaloit to tie 

terminus of each section of inlet manifold, i.e.. 

*CI,J ■ *MI,j * 
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4.5.2 C<dleetor Efftoieney and Collector Outlet Temiwature. tempwattffe of 

eed) eolleetor raust be determined to provide e boundary oonditfcn for the outlet 
manifold heat loss calculations. To accomplish this* the efficient of each collector 
te computed from Equation (2) as follows: 

Th;..!, the collector outlet temperature is computed as follows: 

^CO J “ ^Cl,j ^ 

since 

, . ( 9 ) 

*’1 !,Ag 

4.5.3 Manifold Heat Losses : The equatim for the manifold heat loss was introduced 

previously as 

^LM " ^ V* 

Equation (1) is written with the assumption that the mean value of the manifold flow 
media temperature, tj^, may be computed averaging the temperature at the 
beginning and terminus of the manifold section between adjacent collectors. With 
the appropriate subscripts, may be formulated for both the inlet and outlet 
manifold sections. 

Thus, for the array bank of Figure 3, the mean temperature of the flow media in the 
inlet manifold between collector no. 1 and collector no. 2 may be computed as 

and so on for each inlet manifold lerigth. 
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SSimilar computations may be made for determining the mean temperatiirea, 
in the outlet manifold with the exception that mixing must ocew in t*«e flow 
between the eolleetor outlet ports and the corresponding outlet manifold locatiocu 
Thus* the mixed flow media temperature downstream of the outlet port of c<dleetor 
no. 2 may be congxited as follows 


MOMJ*2 " A 


( 11 ) 


end 80 on for each point at which mixing occurs. The temperature of the flow media 
at the terminus of each section of the outlet manifold may be computed through the 
application of Equation (3) with the appropriate assignment of subeeripts* i.e^ 

^MOM.j • (12) 

tMOB,j - != j =! 5 

4.5.4 Collector Array Bank Effiei«icy : After calculation of the temperature at the 

terminus of the outlet manifold, tj^Q, the collector array bank efficiency and useful 
energy output may be computed from the foUowing equations: 


■ ■ W 

(IS) 

” Qy/dl Ag ^) • 

(14) 


The choice of the gross area, A , to use in these computations depends on the 

o 

effective area for each collector with consideration given to the added area due to 
the external manifold. Thus, the use of Equations 13 and 14 and the appropriate 
effective area permits a comparison to be made with the same array bank, when using 
only the gross area, A^, of the collector alone. 
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5.0 EXAMPLE COMPUTATION 

Two oomm«roially available eollectors wwe selected as examples for tMs analysis. 
The input data was generated from their individual A8HRAB 9S-77 thermal 
performance results and from other assumed parameters listed belowt 


input Paramettf Value fw Each 
Collector 


Inout Parameter 

Definition 

Grumman 

General Electric 

Model No. 

Manufacturer's designu'.ion 

3S2A 

TC 100 

A 

Ccdlector efficiency parameter 

0.730 

0.432 

B 

Collector efficioicy parametor 

0.844 

0.078 

Kar 

Incident angle modifim* 

1.0 

1.0 

^a 

Ambient temperature 

40”f 

40®F 

'mi 

Manifold inlet temperature (40.1,1 50, 220)^F 

(40.1,150,220)^F 

•t 

Total insolation 

«0 

hr.ft* 

sio -S3J 


Gross area of a sinitle collector 

31.80 ft* 

17.41 ft* 

^gEM 

Gross area of a single collector 

35.55 ft* 

19.10 ft* 

including related external manifolds 



Collector fluid flow rate 

400 Ibm/hr 

110 Ibm/ hr 


Mean speci''<% heat of flow 

1.0 

. BTU 

Ibm.^F 

lbm.®F 

n 

Number of collectors in array 
bank 

8 

8 


Area of manifold tube (pipe) 

2.2 ft* 

2.2 ft* 

p 

between collector ports 



R 

Insulation thermal 

(0.75,2.0,7.0) 

(U.75,2.0,7.0) 


resistance of manifold 
covering 

hr.ft*.®F 

BTU 

hr.ft*.®F 

Wtv 


The effective gross area with external manifold was estimated from the 

collector gross area (A^) and the added area due to the manifold and collector 
spacing (Figure 3). 
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The temperatures and total insolation were chosen to produce a typical range of 
thermal performance for the example collectors. The insulation thwmal resistance 

F ^toa 


parameter (R) was varied from essentially no insulation . 

hr ft^ 

relatively high degree of insulation (R » 7.0 ) foe extornal manifolds. 


6.0 KBSULTS 

The manifold heat loss analysis was applied to two separate collectors (Qrummiui 
Model S32A and General Electric Model TClOO) as described in the example 
computation (Section 5.0). The results of the computations are shown in Figures 4 
and 5 for the Grumman and Qmieral Electric collectors^ r^ectiv^y. 


The curves labeled 17 ^ in both Figures 4 and 5 represent the ASHRAB 9S>77 therm U 
performance of the collectors based on the collector's gross area (A^). The curv<is 
labeled are representative of the ASHRA E 93-77 thermcd performance 

modified by the effective gross area (Ag£|^) of the collector including the extern il 
manif<dds. Note that the value of the effective gross area must be careful y 
specified, since it is indicative of the gross "array" area used to produce useful 
collected energy. 


Associated with each curve labeled ^g£|^ are thermal performance curves illustra - 

ing the effects of various pipe insulations. Values of the insulation thermal 

resistance (R) and the associated manifold area (A_) between adjacent collectors are 

P 

input to the computer program. F^res 4 and 5 show the effect on array efficiency 

with R values of 0.75 (essentially a bare pipe), 2.0, and 7.0 * 'Thus, tte 

computer analysis permits optimization of the pipe insulation type and its effective 
thermal resistance. 


Note that the computer printout (see Appendix B) includes all of the inlet and outl< t 
temperatures, manifold heat losses, and efficiencies of each individual collectn*. 
Also, tlie printout includes the array inlet and outlet temperatures, the usefi 1 

t 

energy delivered from the array and the array efficiency. Thus, manifold heM 
losses can be readily evaluated and plotted to indicate the energy lost by the array. 
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CONCLUSIONS 

The computer analysis preser ted in this report permits an evaluation of the effect 
of various manifold insulations on solar thermal coUeotor array efficienev. The 
analysis is simple and produces results quickly. It is believed that the assumptions 
made to develop the analysis aren't serious limitations to the method. The computer 
code can be easily modified. For example, the modification to the code to calculate 
the mean value of the flow media specific heat as a function of temperature may be 
readily accomplished. 

The results indicate significant degradation in collector array thermal performance 
for manifold insulation having thermal resistance values of R - 2.0 hr.ft^.^F/BTU or 
less. Values of R near 7.0 and greater results in little heat loss and performance 
degradation. 

Note that the determination of an effective collector area that includes any 
external manifolding and spacing between a<jyaccnt collectors is critical to adequate 
determination of array performance. The effect of the assumed effective gross 
area for the two example eases in this report show a dramatic effect on array 
efficiency. 
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Manifolds — General Electric Collector, Model TC 100 








APPENDIX A 


LIST OF COLLECTORS IN PHASE I OF THE 
DOE COLLECTOR TEST PROGRAM HAVING 
(INTERNAL) INTEGRAL MANIFOLDS 
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APPENDIX A (eont*d) 

A cursory review of the ooUeotors tested in Phase I of the DOB CoUeetoi* Test Program 
indicate that the list of collectors below have QnternaD integral manifolds. The list 
cannot be cmsidered fully complete or accurate for the following reasons: 

1) Approximately 20 sets of collector drawings and/cr specifications were not 
available for review. 

2) No evaluatimi was made of the flow capacity of the individual collector 
integral manifold, l.e., the number of collectors of a particular modd that 
could be placed in parallel without an external manifold was not determined. 


Ap^oximately 100 sets of collector drawings were reviewed to identify the following 10 
collectors believed to fall in the category of internal manifolds. 

LIST OF COLLECTORS IN PHASE I OF THE DOE COLLECTOR TEST 
PROGRAM HA VINO INTEGRAL (INTERNAL) MANIFOLDS 


NAME 

MODEL 

Daystar Corporation 

21-B, 21-C 

Energy Systems, Inc. 

iiiiD, mis 

Federal Energy Corporation 

p.200 

Gulf Thermal Corporation 

CU30-SL 

International Technology /Solar Corporatimi 

Mark m (40, 
Mark V(30 

Oahu Solar Products 

Solapak 150 

PPG Industries, hic. 

C-524 

Precision Industries, Inc. 

50-11-3MG 

Solar Energy Products, Inc. 

CU-80-WW 

Solar Enterprises, Tnc. 

SLIMLINER 

Solar USA 

28S, 378 

Solaron of North America, Inc. 

SCO-200 

Specialty Mfg. Co. (Insulator) 

SA120-1258S 

Sun Life Solar Products 

SP-100 

Sun Power Systems, Ltd. 

C38B-H 

Sunearth Solar Products Corporation 

3597ADGFB 

Sunworks Division of Enthorpe, Inc. 

N/A 

Wallace Company 

Wallace Company 
Solar Collector 

Wilcox Manufacturing Corporation 

SC8000 


A-2 


APPENDIX B 




SAMPLE PROGRAM INPUT/OUTPUT AND 
COMPOTER LBTINQ FOR COLLECTOR MANIFOLD STUDY 
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APPENDIX B (cont*d) 


SAMPLE PROGRAM INPUT/OUTPUT 
GRUMMAN COLLECTOR, MODEL 332A 
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COifCTOR IRNlfOO) ST(l>v 
l«I£K£PI Eff JCl£NCYs« ;3i 644 

iNLillNl fNi£ NODUlER 6 >lB 26 ,...Se td 

LMI Lm LOW 1066 1066 UHl6 1666 1 066 1 666 1 666 
AClEKi TETfEKRTilf: 46 8 MIFOi) IMIJ lEim(llK£:2^ 6 TGTH Flift>M6 
LOltL W ms f»£6=il 66 COLIEJOR OOSS Mi TH E-ft=35. 55 
CUlECIll: f*6S FUJk («IE=466 86 FJ.U10 bFttlFIC 666 

(f COUiCT05= 8 F»<lFGli) PIPE flIiDKS;* Z 2 mm k W£= L 


ULiiCKfi NO : i 


Ulitl'Ulf NO. : I 


miilUK NO. : i 
UlilCT* NO = 4 


Cai£CTl* NO. = 5 


UlitCIlK NO. : 6 


OliiCI* NO = 1 


auinijk h . : 8 


ICI=*1<).5 ICW25} TflJ=2b,S l(W=i-5 3 TlSii=Z2i 8 

i8ni^«6 dLRinil66 

riuinl=jk66. ri)Olf«^4W Nb=6<24 

TCl=a9 9 ia*^25 2 TRi=a8 S Itti»=225 6 T«re=22l? 

I0FJ=215I.S' feni=<i’4 9 OJUslSri llJI>20i3 
fD0ini=<!bw vjiriF )M) Nb=« ^ 

iu*2w. 6 iopza 1 ini-2iy v 8 iit£= 2 ^. 1 

rWi=2lS8 T6nisi4.8 iuil=lii7 6 (Utt^.3 

mni = 2 m mniriM . Nb=6 224 

TCi=2l?.? TC0=225 6 TW=2a>.7 rWI:224.8 

JBni-218.7 T0PO='.'24.? lfl]=lS7.? aj»=2fc.2 
l|)0l«I=2666 *0Tftl=l666 NG=6 224 

TCi=2l5' 6 TOP224.8 THJ=2lS 6 1/0^4 7 I«06=224 b 

1Bni=2b 6 TBflO=244 . 7 4J!i=l37 8 lUl>26i 1 
dXJTflMb* HwTnCe2666 N0=6 225 

lCJ=2iy 4 JClPZi-i b ini=iH 4 ]fW»=i4. 6 ■,/lje=224 6 

)B«J=21^5 TW=2l'4 6 tUli=l<<7 4 

icoir.i=i266 nb = f > 225 

TLJ=2iy2 ^0^22* S ini^ 2 :tjn=2t’4 5 1W(=224.5 

lWl=<i:y -< )Wl^£i’4 5 {lfll=l87 2 iLtWhi 6 
»0T5J= 866 N0=6 22b 

lLj*2lt * ^'4 1 IB1=216 7 .w^224 4 Inub=o4 4 

IWi=c3\8 ltJlPtl'4 4 ttJ]=18b.8 2Jte2628 

JUtln;* 4t«. f£)0'JlP3266 NO=e 22' 


ccLiCT* fm Tikfonjffta 

T«]J= 226 66 TH< = 224 36 

y(£N= I35t5 »£fla>=6 163 NLW4) 164 


I 



original page is 


B-3 


APPENDIXBJcpnt'd) 


COMPUTER LISTING FOR COLLECTOR MANIFOLD STUDY 


B-4 



1 

m ' 

■» 

m 

ilM 

m 

m 

1N7 

ittS 

1189 

1810 

UU 

1812 

iH2 


. » . ^ I 




'«»•*** T*** ^ 


Di)D67(ii Tnu§).TW(if).iiff!ar- TWfaf>' 

DIHEKKH TCM10..Tlfll(10).TBfV10>,aA(ll).TTO(18) 

««L ll)(IT!D(18),»jrT(ll).IT,IWn^ * 
moTwaD.mim 

awnw - - * « — ~ 

wnF(7,:2K' 

PIWlT(M.t£ W M lUfO? DTTmfnWi^fr ^ ^ ^ 

10Da*)C 

ir«l£l)H)TO» ■ ' 

mW(7.WB7 

nm(' 7W mtc nmcn iid slh «b):t 

l®KX*)ft8 • 

«nE(7,iim ' 

Rimrc ' m in me-m jcnor miF nuFiBrn 
KnF(7,2»f) •- 

FBfrT('Rtfc»a...98WKE5;'D 
BD(S«WT r-"-’ 

«nEf 7 ,lB 5 ) — r— M 

F»FTC' TWIN TIC WIHTTEiigWTaEi'f) • >- 

ENWi*)TI? ■ “ • " — 

wfnrf7,iw) 

FDWRK' TW IN TIC TDfWmif NT TIC '8) " 

WJTEaitfT)" • 

Fwirr(" iwaroLfi mn-'i) 
lODCi tMl '' - - ^ 

i«Tr(7.i00e) 

FCWirT' TWIN DC JWDLNTIIK WlC'f) 

®iKS,*)n 

»JTF(7,1I89) 

F0BfT(' TW IN TIC COLUCTW GROSS NRBl 'fT 
WIX5*)flG 
l«ITEf7,l«0) 

F0WT(' TW IN TIC «0S5 AREN WTH Em.H IWIFOO.'I) 

*ITr(7.10U) 

»F»fn(' TW IN TIC CftlETTOP aON RNTE 't) ' 

• 

F(»fn(' Tw IN w ycriFir ler 'd 
RBIK5*ICP 

*ITEf7,1012) • •• 

F0W)T(' TW IN TIC WN OF TIC WiaLT ‘m-'D 
REA0<S «>AF 

»ITt7.18l<) • • ' 


^ ^ ^ . I 

!■ ^ •■#••• * 


... ^ ^ 


B-5 


;'H/GINAr. PAGE IS 

quality 


1«14 FORHRK' T\fE IN TJ£ fWHFDLf R VRLIE; 't) 
RB»(5-*)R 
IWTF(fc2«l) 

m FORfUTdX, 'CaiiCTW WNIFOD STDPVO 
t«TE(&2W2)flB 

208? FdMRTdX, 'INTERCEFT EFnCIE«V=',F^i 1S(- 
C 'SL0PE='.F5.3) • ■ 



WnE(&20C)»?T 
2083 R*WT(i)C INCIDENT flNGLr roiFIS VfiUB FIR 
C ai8,2fc ..30OE0H5 fll^^^l0F7.7) - 
»I7F«:.20«)m.miT 

20W F(IW?T(lX.'7»IENT7PITOm^ ~ 

c TWJFOP wrr 7ei>EWTiii!r=',F5:i.i5f 

C OOTBL aii^',F5.i) "*■ ^ " 

«ITF(&20e5)«l-flOEIf'' ^ 

2005 FIWf?f(lX-'C(llfrni?G^ ■ 

C 'imro? OIOSS lie BITfl E-?fr^,75. 2) 
wnE(s,20ftb«)OTi;.[? 

m Fornax. 'caifcroR nrss flw Rf?TE=^Ffc Z5n. 

C 'FLUID mine m='.F5. 2) 
WrTE(fc2e07)ND,fP,R 

2007 FDRHfiTdX, 'NUKF (F CaiET05=VI2.Si ' 

C TWIFODPIPEIISKSV-^FSLW 
C 'I«IF01J)RV(IOE--'.F32) 

C IWTIflL C0M)inC»6 

NK1T?II(1>=NC:*«)(IIC - 

TO(1>=(TNIM2 ‘ 

C +(2 *ff^Tfl)/R)A2 *IWTWa)€FW/^^ 

TW(D=TCia> ■ 

Twra:>=aNii+Tcww/2 

DCi100>^Nl 
is)OTHi(;;=(NC-j+i)*fwn: 
?ci(/)=fi«iad)*(2*«)0TwaM^flp/ii) ■ * 

C ♦(2*flP»Tfi)/H)/(2«fr^(J)«(P+flP,^) 
TKI(.D=TCKJ^ 

TBNKJ)=(TMI(HHTCI(J))/2 

eu!ia>=(fF^)*a9fi(j)-TFt) 

100 CONONIC 

«ftT=l 
DO280 

Tro'J'=T[.KJ)+'NG(J;*IT*flG'/'fW^^ 

^')=J*tflOTr 


» ccr;iic 

MHa'^TTPU' 

c 2 *»roTNPa'*CP<flP/«) 

00 W J=2.NT 

TWj'--'rw<H'»iyiTWuM>*Ta>(.n«ff<»Tr) • 

TWi I>=iTWuT)*C2 

-W nriNiF 

0lM«)M.Ni: 

TWtJ^=(TH0H<JWTW!un)-'2 

CfIfTIHlf 


rntTfCR^MT' 

»fw=flifn’aT*flGEW' 
00 W J=l.Nf 


«IT(ci. 3»il )J. TCK J'. T(OvJ5. THI(.n *nHv.T>. TW(J) 
3flte F«Tf Ul>„ 'OlifCTOP NO -■ , I2.3(. 'TCI^',F!i 1. 

c ]x. i. 2 <. '■Tuir'.Fy 1.3X. Wfr ,ni t. 

f ?x. F5D 

t<?i Jtyi' ytd'^ma >, w< jv aah j\ a.H0i j» 

'TPM='. Fy 1, .-5i. 'n«[t r? 1..^. 
f I ? UHtMT.!) 

rflTi*:. M' 'tOTnicn. fi'tW'.ji Nca' 

SRfc -ix injr'. c? i. >. ^ o. 

f NG= .fy 1) 

» ClNT’Hf 

axitny^ 

‘Tiin.-'.r : > 'Jh. = .r :> 

if r»', f. >|M » ji ^ 

. FT 0. . 

C F! ■ :>■ .FT. '' 

GO ’■ • 

iH 


origW''}' 

(if POOR 




\S 

Vi 


H- 


